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The fluctuation- dissipation theorem is used to evaluate the screening factor of nuclear reactions 
due to the electromagnetic fluctuations in a plasma. We show that the commonly used Saltpeter 
factor is obtained if only fluctuations near the plasma eigenfrequency are assumed to be important 
(u) ~ <C T (h = fcs = 1)). By taking into account all the fluctuations, the highly damped ones, 
with u> > u)pe, as well as those with u> < cj pe , we find that nuclear reaction rates are higher than 
those obtained using the Saltpeter factor, for many interesting plasmas. 
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I. INTRODUCTION 

It is known that fusion reactions within a plasma are enhanced. This is due to the fact that ions in a plasma are 
not "naked" , but are surrounded by electrons in the plasma that form a shielding cloud around them. Thus, an ion is 
not repelled, as is a "naked" ion, since the "dressing" of electrons shield the repulsive potential of the ion for distances 
greater than the Debye radius, Re>. The presence of the plasma is taken into account, in what is called the screening 
factor in estimating cross sections of nuclear reaction rates. In the static case, when the velocities of the ions are 
assumed to be very slow compared to those of the screening particles, the screening factor used is the Salpeter factor 
M (valid in the weak screening limit, when Zi^e 2 <C T Rd^e = 1))- Electromagnetic fluctuations are present in 
a plasma, even when in thermal equilibrium. From these fluctuations, we evaluate the screening of the interacting 
particles. Their effect is to enhance the reaction rates, just as a negative potential reduces the Coulomb repulsion. 

We use the formalism of the fluctuation-dissipation theorem to describe the electromagnetic fluctuating fields. 
Using the fluctuation-dissipation theorem, it is possible to study all the fluctuations existent in a plasma. That is, it 
is possible to also include the effect of screening modes that are not eigenmodes (i.e., not at the plasma frequency). 
These modes are highly damped, but although they have a small amplitudes, their overall effect turns out to be great 
because their phase space is very large. In particular, they have non-negligible amplitudes for lu 3> Lu pe , where u> pe is 
the plasma frequency. In estimating the effect of electromagnetic fluctuations on the screening in the limit where only 
frequencies on the order of uj pe are important {lu ~ Lu pe <C T (K = ks = 1)), the usual Salpeter factor is obtained. 

In a non-magnetized plasma, the longitudinal eigenmodes are Langmuir waves with to ~ Lu pe . For usual laboratory 
and space plasmas, this frequency is much less than T and the assumption that w < T is reasonable. However, if we 
are interested in the effect of all the existent modes in the plasma on the screening, it is important not to be limited 
only to modes where lu ~ Lu pe . 

We show, in section II, that the screening can be described in terms of the electromagnetic fluctuations. The change 
in the reaction rates due to the screening is calculated and our screening factor is compared with the Salpeter factor. 
A discussion of our results is given in section III. 

II. REACTION RATES 

The expression for the reaction rate between ions 1 and 2 is § 
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R = I ' dvxf(v{) J dv 2 f(v 2 ) | vi - v 2 | cr(vi,v 2 ) , (1) 

where a is the cross section and f(v±) and f(v 2 ) are the thermal distributions of the particles. The cross section is 
a{vi,v 2 ) — (Sn(E)/E)P(vi,v 2 ), where P(v\, v 2 ) is the penetration factor, E is the kinetic energy of the ions at large 
separations in the center of mass reference frame and Sn(E) is a purely nuclear factor that varies slowly with E. P(v±, 
v 2 ) is proportional to exp[(— 2(2fi) 1 ^ 2 )/h) f r ° [U(r) — Uq — E^^dr], where fx is the reduced mass. The integration is 
performed from the nuclear radius a to the turning point radius ro. U(r) is the bare Coulomb interaction ZiZ 2 e 2 /r 
and Uq, the potential energy due to the screening. In the case of weak screening, Uq is piratically independent of r 
in the integration interval. The reaction rate R is then found to be increased by a factor F = exp(Uo/T) 0. In the 
static case, the weak screening limit, the enhancement factor is exp{Z\Z 2 e 2 / (T Rd)), the Salpeter factor. 

In a plasma, electromagnetic fluctuations are always present, acting as a background for the interacting ions. Let us 
consider the effect of the fluctuating electric fields. <f>, the electric potential acting in the region from ro to a, creates 
a potential energy Ui — e(Z\ + Z 2 )<j). The enhancement factor is, then, F — exp(Uo/T) = exp(U\ — Uao/T), where 
Uoo is the potential energy of the ions when they are at infinity, Uoo = Zie0i jOO + Z 2 e<fi 2 ,oo. Averaging F in time and 
using a Taylor expansion, we obtain terms which are proportional to (4> 2 ), ((f) 2 ^) and ((f> 2 OQ ), as well as cross terms. 
Because the fields 0, 4>i.oo an d 02, oo are random with respect to one another, the only terms that remain to first order 
are (<fi 2 ), (<j>\ ^ and (cf> 2 QO ). Assuming that they are equal to one another, 

F = 1 + ZiZ 2 e 2 (0 2 )/3 2 + ... = exp (ZiZ 2 e 2 (0 2 )/? 2 ) , (2) 

where /3=1/T. 

4> is assumed to be changing slowly in the integrand of the penetration factor in the region from r to a. Therefore, 
we consider fields with frequencies to < uj ma x = 2irro/v and wavenumbers k < k max — 2it/rQ, where v is the velocity 
of the reduced mass in the center of mass system. We find that (<fi 2 ) is relatively insenstive to the exact upper limits 
of the integrations, uj max and k max . 



b 2 ) is related to the fluctuations of the longitudinal electric field in a plasma: (<j> 2 ) k = (E 2 ) k /k 2 . Therefore, 

/ " n ' duo^-^ . (3) 
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^From the fluctuation-dissipation theorem [y|JJ] , the expression for the intensity of the fluctuations of the longitudinal 
electric field is 
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where Si is the longitudinal dielectric permittivity of the plasma. We note that this is an exact quantum mechanical 
relation. Given a dielectric permittivity that describes the plasma, it is possible, from Eq. (4), to obtain the 
fluctuations in the longitudinal field due to the eigenmodes (i.e., the electrostatic Langmuir waves), as well as to all 
the other fluctuations existent in the plasma. 

We show below that the overall effect of fluctuations with ui > cu pe is to substantially increase the reaction rates. 
Although the amplitudes of these fluctuations are small, they are not negligible. Due to the very large range of 
their frequencies, their overall effect is to increase the reaction rates appreciably. Because the reaction rates are 
exponentially dependent on the screening potential, its increase is appreciable for reactions with large values of Z\Z 2 
such as p — 7 Be and p — 14 A. In the limit where only fluctuations with frequencies on the order of the eigenfrequency 
Lu pe (<§; T) are important, Eq. (4) can be easily intergrated over frequency with the aid of the Kramers-Kronig 
relations, obtaining 

<■■>•- 1 C 1 - sob>) ■ < 5 > 

This is the expression which is generally found in text books &£§. For an electron plasma in the collisionless case, 

K L 
k- 
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the longitudinal dielectric permittivity is £/(o;,k) = f + -£■ ( 1 — <fi(z) + i^/irze~ z ), where z = ^3/2(uj/kvT) and 
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2 z e z J e x dx. With this dielectric permittivity, (<fi 2 ) = T/Rd [Eq.(3)] and the enhancement factor F, in 



Eq.(2), is the Salpeter enhancement factor, 
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We use the model described in detail in Opher & Opher which includes thermal and collisional effects, to 

obtain the longitudinal dielectric permittivity. It uses the Vlasov equation in first order with the BGK (Bhatnagar- 
Gross-Krook) collision term, which is a model equation of the Boltzmann collision term @] . The inclusion of collisions 
changes the results very little and the effect on ej(o;, k) for the plasmas studied, is almost negligible. For example, 
in the case analysed in Opher and Opher M, where T = 10 5 K and n = 10 10 cm~ 3 , the difference in the correlation 
energy, with or without collisions, is less than 10~ 6 . 

The dielectric permittivity that we therefore use, is almost identical to the collisionless dielectric permittivity given 
after Eq. (5) in the text (the only difference is the inclusion of other shielding species besides the electron). From 
this model, the longitudinal dielectric permittivity for an isotropic plasma is found to be 



(u+iri) y f 0)+ 



2 1 + ^!L>Z 



£/(cj,k) = 1+ > : — - , (7) 

^ V2kv a V ^kv a 

where a is the label of the species of the particles, v a is the thermal velocity of the species and Z(z), the Fried and 
Conte function. By using Eq. (@) in Eq. (§) and Eq. (@), we obtain (</> 2 ). The enhancement factor is obtained from 
Eq. (2). We note that using a relation which is similar to Eq. (2) for the transverse field (B 2 ) kyUI in conjunction with 
the transverse dielectric permittivity instead of Eq. (7), the correct black body spectrum is obtained [5,6]. 

We estimate the enhancement factor F for the reaction rates p — p, 3 He — 3 He and p — 14 N for an electron-proton 
plasma at temperatures T = 10 7 — 10 8 K and densities ranging from 10 23 — 10 28 cm -3 , choosing the maximum density 
for each temperature so that the plasma parameter is less than g = l/n\ 3 D ~ 0.8. When g approaches such high 
values, non-linear effects begin to be important. (Note that Ad has an extra factor V2 due to the inclusion of protons.) 

The Gamow energy, the most effective energy for thermonuclear reactions, is Eq = \.22{{ZiZ2) 2 ATq) 1 ^ 3 keV, where 
A is the reduced atomic weight and Tq = T/10 6 . Therefore, for each of the temperatures that we choose, the Gamow 
energy changes, with a consequent change in the limits of the integration, k max and uj m ax- As a particular case, we 
assume the center of mass energy to be zero. (F — Fsal)/Fsal(%) vs g, where F$al is the Salpeter enhancement 
factor and g is the plasma parameter, is plotted for T = 10 7 K in Fig. 1. In Fig. 2, we plot (F — Fsal)/Fsal{%) vs 
g for T = 10 s K. We also estimate the enhancement factor F for the reaction D — D in an electron-deuteron plasma 
for temperatures T = 10 7 — 10 8 K and densities 10 23 — 10 28 cm~ 3 . Again, the maximum density for each temperature 
is chosen to be g ~ 0.8. (F — Fsal)/Fsal(%) vs g is shown in Fig. 3. We see that the enhancement factor F is 
larger than the Salpeter factor. 

In order to see how the fluctuating potential is affected by these modes, we plot the potential (0 2 ) w k vs ui/uj pe for 
the reaction p — 14 N, with T = 10 7 K and n = 10 24 cm -3 (in this case, g = 0.3) in Fig. 4. The solid curve is plotted 
for k — 100 fee, the dashed curve for k = 10 fco and the dotted curve for k — ko- We normalize the curves and plot 
(</> 2 ) w k divided by ((j> 2 ) k w=0 - When (<fr ) u k is integrated for all wavenumbers k, but only for frequencies up to ~ aj pe 
(<C T), we obtain the Salpeter potential. We see that for k — ko, the potential ((f) 2 )^^ has non-negligible values only 
at low frequencies. However, for k = 10 ko and k = 100 ko, the potential cj) k ^ spreads out to high frequencies. For 
these high values of k, instead of dropping abruptly at u> ~ uj pe , the curves decrease slowly with frequency. There 
is thus a substantial contribution for u> > to pe . Thus for large wavenumbers, (i.e., k > ko) fluctuations make a 
substantial contribution at frequencies higher than the plasma frequency. For example, for k = 10 fco and u) — 20 oj pe , 
(^ 2 )u),k/((f> 2 )uj=Q,k is not zero, but 0.002. Similarly, for K = 100 k D and uj = 20 uj pe , (0 2 ) w ,fc/(0 2 )w=o,fc is 0.01. 



III. CONCLUSIONS AND DISCUSSION 



By using the fluctuation- dissipation theorem to estimate the screening of the reacting ions due to the electromagnetic 
field fluctuations in plasma, we calculate the enhancement factor for the reaction rates. The enhancement is found 
to be greater than that for the adiabatic case (u) -C T), where the Salpeter factor is valid. It is important to note 
that this effect is not due to a specific choice of the dielectric permittivity. The model that we use includes the 
standard dielectric permitivity and gives almost the same results as does a collisionless thermal description. A higher 
enhancement factor is obtained because we include the high frequencies and do not assume that lu < uj pe <§; T. Using 
this method in a previous study to derive the transverse magnetic fluctuations, the standard black body spectrum 
was obtained. For the plasmas studied here, electron-proton plasmas with temperatures 10 7 — 10 8 K and densities 
10 23 — 10 28 cm -3 , the reaction rates for p — p, 3 He — 3 He and p — 14 N are increased up to 8 — 13% in the regime 
g < 1. The reaction D — D for an electron-deuteron plasma is increased by 1 — 2%. We find that these fluctuations 
have small, but not negligible amplitudes, compared to those at low frequencies. They exist, however, in a much 
larger range of phase space, up to very high frequencies. Their overall effect is to increase the screening potential. 
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Because the reaction rates are exponentially dependent on the screening potential, this increase is especially strong 
for reactions with relatively large Z\Z<i, such as p — 7 Be and p — 14 N. When summing the fluctuations only up to 
the plasma eigenfrequency ui ~ u> pe , the Salpeter enhancement factor is obtained. We show here that the fluctuating 
potential has a large contribution from highly damped modes with frequencies u> > Lo pe . 

There has been some controversy about the existence of a screening factor, which is different from the Salpeter 
factor, due to dynamic screening ||. The effect described here is due to quasi-fluctuations that are heavily damped 
and that always exist in a plasma . They alter the Salpeter screening factor, although not necessarily in the same 
way as does dynamic screening. 

This work shows that by taking into account all the modes, including those which are heavily damped, with 
frequencies higher than the plasma frequency, the screening factor of the nuclear reaction rates is altered and is larger 
than the Salpeter factor. The consequences of this work are left for future studies. 
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Department of Energy for support (grant De-FG03-93ER54224). 
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FIG. 1. The expression (F — Fsal) / Fsal (%) vs g for an electron-proton plasma at a temperature T = 10 7 K. F is the 
enhancing factor, which is compared to the Salpeter enhancing factor Fsal- g is the plasma parameter l/n\ A D . The solid curve 
is the p — p reaction, the dashed curve is the 3 He — 3 He reaction, and the dotted curve, the p — 14 N reaction. 



FIG. 2. The expression (F — Fsal) / Fsal (%) vs g for an electron-proton plasma at a temperature T = 10 s K. F is the 
enhancing factor, which is compared to the Salpeter enhancing factor Fsal- g is the plasma parameter l/n\ 3 D . The solid curve 
is the p — p reaction, the dashed curve is the 3 He — 3 He reaction, and the dotted curve, the p — 14 N reaction. 

FIG. 3. The expression (F — Fsal)/ Fsal (%) vs g for an electron-deuteron plasma. F is the enhancing factor,which is 
compared to the Salpeter enhancing factor Fsal - g is the plasma parameter ljn\ 3 D . The solid curve is the reaction D — D for 
T = 10 7 K and the dashed curve, the reaction D - D for T = 10 8 K. 

FIG. 4. The potential ((j> 2 ) w k/( ( t )2 )u:=o k vs W^pe f° r the reaction p — 14 A in an electron-proton plasma at T = 10 7 K and 
n — 10 24 cm -3 . The solid curve is for k = 100 kn, the dashed curve is for k = 10 kn, and the dotted curve, for k = ko- 
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